Lecture 6 

The 3D Gasket 

Warning: this is a draft copy. It has not been passed any revision 



How to construct the 3D Gasket 



We start with a tetrahedron 



We thin find and connect the six midpoints on the edges of 
the tetrahedron. This gives four smaller tetrahedrons, one 
for each original vertex and a cutout area 



Recursively continue subdivide each resulting tetrahedron 
for a specified numbers and finally draw the resulting small 
ones. 
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3D Gasket program 



We need a function to draw the 
resulting tetrahedron triangles 



We need a function that uses the 
above function to draw a 
tetrahedron 



]void triangle (GLfloat *va., GLfloat *vb, GLfloat *vc) 
{ 

glVertex3lV( va) ; 
glVertex3f ^(v/b) ; 
glVertex3f v( vc) ; 

> 

]void tetra (GLfloat *a., GLfloat *b, GLfloat *c, GLfloat *d) 
{ 

glColor3fv( colors [8] ) i 
triangle(a, b, c); 
glColor3fv(colors[l] ) ;| 
triangle(a ± d); 
glColor3fv(colors[2] ) ; 
triangle(a, d, b); 
glColor3fv(colors[3] ) ; 
triangle (bj d^ c) - 
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3D Gasket program: division of the tetrahedron 



]void divide_tetra(GLfloat *3 J GLfloat *b, GLfloat *c s GLfloat *d_p int m) 

{ 

GLfloat mid[6][3]; 
int j; 
if (m>&) 

{ 



else( tetra(a ± b,c ± d) ) ; /* draw tetrahedron at end of recursion */ 



for(j=&; j<3; mid[&] [ j ]=(a[ j ]+b[ j ] )/2; 

for(j=*j j<3; mid[l][j] = (a[j]+c[j])/2j 

for(j=0j j<3; mid[2][j]=Ca[j]+d[j])/2j 

for(j=*j j<3; mid[3] [j ]=(b[ j ]+c[ j ] )/2j 

for(j=0j j<3; mid[4][j]-(c[j]+d[j])/2j 

for(j=0; j<3j mid[5][j]=(b[j]+d[j])/2j 



divide_tetra(aj mid[0]., inid[l] f mid[2]j m-l)j 

divide_tetra(mid[9] t b t midfS], mid[5] t m-1); 

divide_tetra(mid[l] , mid[3], c, mid[4], 

divide_tetra(mid[2] t mid[4], d, mid[5]j m-l)j 




/* create 4 tetrahedrons by subdivision */ 




:■ 
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3D Gasket program: integrating the parts 



GLfloat v[4][3]={{9.e ± {&.&, &.942S&9, -&. 33333}, 

{-0.S16497, 471405., -0-333333}., {0.S1S497, -0.4714&5, -0-333333}}j 

GLfloat colors [4] [3] = &-0}., {0-0, 1-0, &.0}j 



void display() 
{ 

glClear<GL_COLOR_BUFFER_BIT | G L_DEPTH_BUFFER_BIT) j 
glBe gin (GL_TRI AN GLES) j 

divide_tetra(v[&], v[l], v[2] t v[3] t n)j 

glEnd(); 

glFlush()j 
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Hidden surface removal 



glClear (GL_CQL0R_BUFFER_B1T I GL_DEPTH_BUFFER_BIT) ; 

GL_DEPTH_BUFFER_BIT makes OpenGL clear the z 
buffer in the frame buffer, when clearing the graphics 

glutlnitDisplayModeCGLUT.SINGLE I GUTT.RGB I CLtTT.DEPTH) | 

GL_DEPTH makes OpenGL to enable the z or depth 
buffer (a layer in the frame buffer) 



int main(int argc., char **argv) 



glEnable (GL_D£PTH_TEST) 



n=atoi ( argv[ 1 ] ) j /* or enter number of subdivision steps here */ 
glutlnit(&argc, argvji 

glutlnitDispl ayMo d e (GLUT_SI MGLE | GLUT_RGB | GLUT_DEPTH ) ; 

glutInitWindowSize(5&& J 509); 

glutCreateWindow( "3D Gasket" ); 

glutReshapeFunc(myReshape) ; 

glutDisplayFunc(display) ; 

glEnable(GL_DEPTH_TE5T); 

glClearColor (1.©, I.e., I.e., 1-9); 

glutMainLoop( ) ; 



GL_DEPTH_TEST Make OpenGL enable the 
hidden surface removal algorithm that uses 
depth (z) to remove hidden surfaces 
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Drawing Implicit and explicit functions 

Implicit function: g(x,y) = 0 Explicit function: Y = h(x) 

Emplicit example (circle centered at the origin): x 2 + y 2 — 1= 0 
Emplicit example (Ovals of Cassini): (x 2 + y 2 + a 2 ) — 4a 2 x 2 — b 4 = 0 

Explicit (Sine function): y = sin(x) 



Explicit functions are already separable. Then we can substitute in the function 
with the input and draw 

In the case of Implicit function: The function may generate no, one or many 
curves depending on the parameter and the equation (example 2, for example, 
may generate no, one or two curves depending on the values of a and b 
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Explicit function: General form 

General: Z = /(x,y) Contour: c = f(x,y) or f(x,y) — c = 0 



In this general form, instead of having no, one or many curve for the case z=0, 
we will have all these possibilities for every possible value of z. The curves 
generated for a given value of z=c are called the counters for z=c. 



To find the contour for c, two ways could be taken: 

• Solving the implicit equation: f(x, y) — c = 0, for the x , y values (difficult 
or impossible, no systematic way to go on) 

• Sampling the region of interest in both x and y, study the function values at 
these samples to see if the contour pass by them or not. If it passes draw it, 
otherwise don't draw, (there is a systematic way to go on) 

The second way is an imitation to sampling an object using devices such as 
laser or satellite. Note also that the second way only approximate the contours 
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Approximating a contour: Marching squares 



We start by a grid of x and y (points where samples will be tested). For each 
cell of the grid, we see if the contour pass by the cell or not. If it passes, how it 
passes (approximation of the intersection point(s) of the contour with the cell 
edges. Finally, we connect these points of intersection using line segment to 
get an approximation for the contour. The following is an example for a cell in 
the grid and the contour: f(x , y)=c 
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For these conditions, we 
have the possibilities: 




We approximate or interpok..^ 
based on the values of f(x , y) 
at the vertices assuming single 
intersection 



r 
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Intersection possibilities 



There are 16 (2 4 ) way to color the vertices in black (less than c ) and white 
(greater than c). From these 16 there are only 4 unique cases. Each of the 
other cases is symmetric with one of these 4 
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What we need now is how to draw a line 
segment inside each cell (if one exist) to 
approximate the contour. The line segment 
could be unambiguously approximated 
except in the last case which suggests two 
ways to put the two line segments (a, b). 
Which pair we choose: Random choice, 
Always use one of them, or subdivide the 
cell to remove ambiguity 



O 
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